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Abstract: We study the problem of searching for cosmic string signal patterns in the 
present high resolution and high sensitivity observations of the Cosmic Microwave Back- 
ground (CMB). This article discusses a technique capable of recognizing Kaiser-Stebbins 
effect signatures in total intensity anisotropy maps from isolated strings. We derive the sta- 
tistical distributions of null detections from purely Gaussian fluctuations and instrumental 
performances of the operating satellites, and show that the biggest factor that produces 
confusion is represented by the acoustic oscillation features of the scale comparable to the 
size of horizon at recombination. Simulations show that the distribution of null detec- 
tions converges to a 'x^ distribution, with detectability threshold at 99% confidence level 
corresponding to a string induced step signal with an amplitude of about 100 /uK which 
corresponds to a limit of roughly G/i ~ 1.5 x 10~^. We implement simulations for deriving 
the statistics of spurious detections caused by extra-Galactic and Galactic foregrounds. 
For diffuse Galactic foregrounds, which represents the dominant source of contamination, 
we construct sky masks outlining the available region of the sky where the Galactic con- 
fusion is sub-dominant, specializing our analysis to the case represented by the frequency 
coverage and nominal sensitivity and resolution of the Planck experiment. As for other 
CMB measurements, the maximum available area, corresponding to 7%, is reached where 
the foreground emission is expected to be minimum, in the 70-100 GHz interval. 
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1. Introduction 

Current theories of particle physics predict spontaneous symmetry breaking in the early 
Universe, with the consequent creation of coherent structures in the spatial distribution of 
the fields involved in the process; these structures, known as topological defects, represent 
genuine tracers and unique remnants of those processes, occurring at energies which are 
not achievable with the current terrestrial particle colliders, see Q and references therein. 
Among the most known and studied relics of this kind, cosmic strings correspond to an 
uni-dimensional region in which a (scalar) field is trapped by true minima of its potential, 
storing energy in that region. The energy density is parametrized by the dimensionless 
quantity G/i where G is the Newton's constant, and /i the string tension, see Q for review. 
Years ago, while it was still unclear if strings could dominate the structure formation pro- 
cess, their signal was searched in the main statistical indicators of structure formation 
itself, like the power spectra of anisotropies in the Cosmic Microwave Background (CMB) 
or Large Scale Structure (LSS). The conclusion of these studies was that strings should 
contribute to cosmic structure formation by no more than a few percent ^ ^ |3^, ^ . 

Indeed, the evidence for coherent acoustic oscillations in the CMB and LSS power spectra, 
dominated by density fluctuations, favored the scenario of isentropic or adiabatic Gaussian 
fluctuations like predicted in the simplest inflationary scenarios, and indicated that cosmic 
strings, if any, played a minor role in providing the initial conditions, in a statistical sense, 
to the Universe we see today. Early in their study, it was evident that strings produced 
perturbations that were not coherent as the observed acoustic oscillations and that they 
tend to produce equal power in scalar (density), vector (vorticity) and tensor (gravitational 
waves) perturbation modes. Recently, |Q proposed to look at sub-dominant, but most 
interesting components of the power spectrum of CMB anisotropies, namely the B modes 



activated by cosmological gravitational waves, to constrain cosmic strings. 
The hypothesis is then that these objects are washed away by inflation itself motivating 
the direct search of the signal from isolated and rare cosmic strings. Among these studies, 
those on CMB are based on the well known Kaiser-Stebbins effect |1S] resulting in a step 



like feature in the CMB temperature caused by a string which is orthogonal and moving 
relativistically with respect to the line of sight. These studies produce upper limits on the 
abundance of string in our own observable Universe of course, but more quantitatively on 
the string density parameter G/i |16, 0, |2^, 29, 20, |2|, On the other hand, studies 



in the optical band look for strong lensing events, in which the strings bend the light from 
a distant Galaxy, causing a mirror image of that |Q . 

In this paper, we speciflcally search for signals imprinted on CMB temperature observa- 
tions by long cosmic strings of cosmological scale, which was predicted by theory and also 
supported by simulations^. The long strings between the last scattering surface and us 
would be back-lighted by CMB photons and leave signatures in CMB temperature map via 
the Kaiser-Stebbins effect. While the works mentioned above are dealing with the effects of 
cosmic string network on statistical observables, this work aims at the direct search of sig- 
nals produced by isolated cosmic strings left on CMB temperature anisotropy maps, having 
a coherence length corresponding to at least the angle subtended by an Hubble horizon 
at decoupling. The Planck survey^ of CMB temperature and polarization anisotropics is 



ongoing |24|, promising a leap forward in many research fields, including constraints on 
cosmic strings abundance and tensions. The unprecedented frequency coverage, sensitivity 
and angular resolution requires a dedicated discussion targeted to the applicability and 
expectations from Planck direct searches of cosmic strings. In this paper we work out two 
missing issues concerning the direct search of KS effects of cosmic strings. The first one is 
represented by the statistics of null detections from pure Gaussian CMB anisotropics, with 
nominal noise and angular resolution of the operating satellites, which, needs to be quan- 
tified prior to data analysis. Consequently, as the second step, one can use the available 
models and data of Galactic and extra-Galactic foreground signal simulations at Planck 
frequencies, to select the sky area where false detections provide a negligible effect with 
respect to the main contaminant represented by the background, Gaussian and uniformly 
distributed, CMB anisotropics. In this paper we discuss and specialize the technique dis- 
cussed in [0] to deal with these issues. In Section 2 we review in detail the effects of cosmic 
strings on the CMB temperature anisotropics and outline the string search algorithm. In 
Section 3, we apply the technique to simulated Planck data sets, quantifying the statistics 
of null detection based on the nominal performances of the satellite, as well as selecting 
the sky region which is available for searches based on the current knowledge of Galactic 
and extra-Galactic foreground emission. Concluding remarks are in Section 4. 

2. Cosmic String effect on CMB and detection algorithm 

Consider a cosmic string with mass per unit length /i, velocity /3 and direction s both of 



^http:/ /www. damtp.cam.ac.uk/research/gr/public/cs_evol. html 
^ www.rssd.esa.int / planck 




Figure 1: Decomposition of a patch of sky map which contains a discrete temperature step. The 
gray spots represent masked point sources. The discrete temperature step is slanted 15° from 
the polar axis. The instrumental noise was generated with power <jn = and the height of 

temperature step is lOO/xK. |^J.J 

which are perpendicular to the line of sight, back-lighted by a uniform black body radiation 
background of temperature T. Due to the angular deflection 5 = SvrG/u in the conic space- 
time around the string, there is Doppler shift on photons passing around it on different 



sides, causing a temperature step across the string [19, 03] with magnitude given by 



signature through {2A_) is simply the angle subtended by the CMB horizon at decoupling, 
as the maximum extent of the conic space-time formed around the string would be limited 
within the horizon at decoupling. The largest causally connected patch possible would 
be optimal in this pattern search, because a gradient or step-like pattern generated by 
white noise would be suppressed by ~ 1 / ^ iVpixei factor, and our choice the size of a test 
patch as the horizon size at the decoupling can be justified. The anisotropy in the CMB 
temperature sky map is caused by several processes including intrinsic perturbations in 
the photon fluid, metric fluctuations, possibly topological defects like cosmic strings. CMB 
photons pick up secondary anisotropies when traveling across forming structures on the 
line of sight toward us. Foreground anisotropies are caused by the emission from extra- 
Galactic point sources and our own Galaxy. Finally, the instrument adds instrumental 
noise and various systematics. By considering intrinsic CMB anisotropies, noise and the 
string contribution, one has 

2^pixcl = TCMB + Tnoise + ^step • (2.3) 



o 
o 



^ = 8^G^7/3 , (2.1) 



o 
I— ^ 
o 



where 7 = l/\/l — In a matter dominated universe the projected angular length of 
string in the redshift interval [zi , Z2] scales as ^/zl - The standard cosmological model 
{^Ide ~ 0.7, 0,m ~ 0.3) places the decoupling of the CMB, also known as last scattering, 
at an epoch which corresponds to a redshift zis ~ 1100. The apparent angular size of the 
CMB horizon at last scattering is therefore given by O 

1 /iooo\^/^ O 

R = ^ = 1.8° ^ ~ 1.7° . (2.2) ^ 

Thus, a convenient choice for the angular size of patches to be examined for cosmic string 
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We employ five parameters which characterize a circular patch of string embedded sky 



map, To, (Jg, A, p and a [18|. They represent the background uniform CMB temperature, 
the sum of the CMB anisotropy and noise variances, the amplitude of the discrete temper- 
ature step, the ratio of the blue-shifted pixels with respect to the total, and the orientation 
of the step, respectively. To recover these parameters from a given CMB region to be 
examined, where the individual components are superimposed, we define five observables 
which can be expressed in terms of (To, ug, A,p, a): 
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where the integrations are done over the examined sky area, chosen circular for simplicity, 
with radius R corresponding to the angular size of the horizon at decoupling, as we men- 
tioned above. When a discrete temperature step is present and all the other components 
are suppressed, the observables defined in ( |2.4| ) - can be analytically computed and 
expressed in terms the five step parameters: 
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where the discrete temperature step can be written as 



f step 



' (x - Xc) - - 



(2.12) 



and Xc {—R < Xc < R) is the x-coordinate of the step with x-axis running normal to 
the temperature step, and 9 {x — Xc) represents a step function. The relation between p 
(red-shifted area/total area) and Xc (coordinate of discontinuity) is given by 



1 

P=2 



1 



1 - - - sin"^ ( — 

R^ TT \R 



(2.13) 



The expressions for the observables (U) - (|1|) can be used to find the step parameters 



by solving for them, recovering the input ones. Now, taking into account the presence of 
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the non-string CMB fluctuations and noise, the observables can be calculated to the first 
order in the step parameters: 

- c^cMB + ^^ + A2p(l-p) 



a 
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(2.14) 
(2.15) 
(2.16) 
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(2.17) 



At this point, we solve for the temperature step parameters (a (orientation), Tq (uniform 
background temperature), Xc (location of temperature step), A (height of temperature 
step), aQ = o"cMB + '^n) ™ terms of observables: 

'D^+iDy 
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(2.18) 
(2.19) 

(2.20) 
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2 I 2 
<^CMB + 



a 



(2.21) 

A^p{l-p) . (2.22) 

Note that we use p and Xc interchangeably since they are related by ( |2.13| ). Expressions 
(2.18) - ( 2.22| ) are the key formulas in the algorithm used to compute the step parameters 
in any give horizon sized circular patch of the sky. Due to the presence of CMB and noise 
fluctuations (and foreground emissions, to be considered next), the results will have statis- 
tical fluctuations. The dispersions of the outputs are the estimate for the error associated 
to the quantities measured by the algorithm discussed here. 
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3. String detection statistics and foreground effects for satellite measure- 
ments 

In this section, we apply the search algorithm outlined in the previous section in order 
to examine two important aspects of string search for satellite CMB measurements. The 
latter are characterized by a wide sky area and frequency coverage, mapping in particular 
the Galactic foreground contamination to high accuracy in bands where it is dominant, 
as well as high angular resolution and sensitivity. When looking for string signatures, one 
has to characterize the statistics of null detections caused by a pure Gaussian spectrum of 
CMB anisotropics, in order to assess the probability that a given event is actually caused 
by a string and not by an ordinary CMB fluctuations; moreover, the spurious detections 
caused by foreground emissions have to be quantified and controlled. In this Section we 
address both aspects. 
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67 
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29 


41 
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Table 1: Comparison of cumulative occurrence levels of confusion outcomes for large and small 
scale CMB anisotropy power, with varying noise amplitude, map resolution. 



3.1 Null detections from Gaussian CMB anisotropies 

We use the publicly available simulation codes CAMB'^ and HEALPix^ |12] to generate 
simulation maps with a consensus A Cold Dark Matter (ACDM) cosmology input param- 
eters = 0.046, QcDM = 0.233, Qa = 0.721, h = 0.701) (see |2| for updated results 
on cosmological parameters). We run the detection algorithm described in the previous 
Section on these maps, and evaluate the statistics of null detections, adopting the nominal 
performances of the currently operating WMAP and PLANCK satellites. Indeed, we do 
not expect significant difi^erences in the two cases, as the signal we are looking for is mainly 
concentrated on the degree scale, where both instruments are substantially limited by cos- 
mic variance. The results are shown in Figure |2|, solid line with colored integral. On the 
left, the angular resolution and nominal sensitivity were set corresponding to the WMAP 
W band with central frequency at 94 GHz, FWHM=13.2 arc minutes, A'^^ide = 2^, white 
noise noise rms = gm = 35//K on 0.3° x 0.3°. On the right, the angular resolution and sen- 
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sitivities were chosen correspondingly to the PLANCK 100 GHz channel on-board the High 
Frequency Instrument, with FWHM=10 arcmin, A^side = 2^^, noise rms = ajsf = 2.5/iK on 
squared pixels with side equal to the FWHM. The first thing to note is that the distribu- 
tions of the null detected A^^^. are similar as we anticipated, and approximately follow a 
X^-distribution, due to the randomness of Gaussian peaks which determine the spurious 
signal. The color filled curves in Figure ^ show the different confidence level that a given 
event in a real measurement is not due to spurious detection from Gaussian anisotropics. 
For WMAP , 68% of occurrences falls in < A < 77/iK, 95% in A < llO^uK and 99% 
in A < 147^K, respectively. For PLANCK , we have 68% in < A < 75/iK, 95% in 
A < 117/iK and 99% at A < 143/iK. The dotted lines also show the cases in which the sky 
signal was added with a step with Ain = lOO^K in each area examined, clearly indicating 
that the search algorithm is working, as the distribution peaks at the chosen value of Ajn. 
The fact that WMAP and PLANCK show nearly same results indicates that the effects of 
map resolution or instrumental noise level applied here on the result are nearly negligible 
compared to other factors, which we will discuss further in the following. 
In conclusion of this Section, we show that the main contaminant for cosmic string search 
with specifications corresponding to the operating satellites are not given by the noise but 
by the underlying Gaussian fiuctuations, which are dominated by acoustic oscillations on 
the degree scale which we are focusing on. In order to contrast the effects of super-horizon 
scale fluctuations and sub-horizon ones, we compute the statistics of null detections by con- 
sidering simulated maps with anisotropy power on the super-degree angular scales (/ < 140) 
and sub-degree ones (/ > 140), separately. Notice that the rms fluctuation power is about 
the same for the two sets: 

1=2 i>140 

In Table |, one sees the variation of the confidence levels for the large and small scale power 
as a function of map resolution, noise level. At the first glance, it is clear that, while noise 
and map resolution do not induce significant differences, the most significant improvement 
is obtained when super-degree and degree scale power is taken out. 

3.2 Spurious string signals from foregrounds 

Foregrounds are considered as one of the ultimate limitation of CMB measurements, and 
represent a major source of contamination, in particular for cosmic string searches. Multi- 



frequency data analysis techniques are being studied for PLANCK [22|, and applied suc- 
cessfully to the WMAP data in order to separate them from the main CMB emission. 
Here we take a conservative approach, by exploiting the current foreground observations 
outside and inside the microwave band, performing an all sky analysis looking for regions 
where they cause spurious string detections, focusing on the PLANCK frequency channels 
going from 30 to 353 GHz. 

Let us consider extra-Galactic point sources first. The usual technique to identify and 
remove them is to isolate 5a events above the CMB and noise standard deviations, and 
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30 GHz 


69 


109 
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70 GHz 
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100 GHz 


77 


119 


147 
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77 
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77 


119 


147 
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77 


119 


147 



Table 2: Cumulative occurrences of string signal detection from the simulated full sky CMB maps 
of 7 PLANCK frequencies. For each frequency band, smoothing and instrumental noises are applied 
accordingly as specified in pTt - 



removing them with optimal filters [^] . Residual sources behave as a correlated noise due 
to beam smoothing which remains in the maps. We simulated randomly scattered point 
sources according to the existing modeling of the radio and infrared populations H, 11 1. 
After removing the 5a signals, causing typically a few holes of beam size in the patches 
we consider, we run our cosmic string search algorithm to look for spurious detections. 
By comparing the recovered spurious A, we find that no significant disturbance from un- 
resolved point sources is observed, due to their sub-dominance with respect to the CMB 
acoustic oscillations, which largely prevails as a contaminant. 

The situation is markedly different for the diffuse foregrounds coming from our own Galaxy. 
We construct all sky maps including the effect of three of the main diffuse Galactic emis- 
sion mechanisms: the synchrotron emission is caused by electrons spiraling around the 



lines of the Galactic magnetic field, using the data by [y^, treated by [10| with spectral 
index obtained from WMAP analysis [|l|, see for more detail^; the dust total intensity is 
based on the analysis of IRAS and DIRBE data by 0] , implementing model 8 of frequency 
scaling including spatial variations of dust frequency scaling; the free-free emission from 
Brehmstraahlung of electrons hitting ions in the Galaxy, traced by Ha emission, has also 
been included^. 

We simulate sky maps with these emissions included at 30, 44, 70, 100, 143, 217 and 353 
GHz, and apply the cosmic string detection algorithm. As expected, we do find regions 
with high rate of contamination, concentrated close to the Galactic plane, at all frequen- 
cies. We show the results in Figure ^, where we show all together the regions yielding 
spurious cosmic string signals with various confidence levels based on the distribution of 
null detection from CMB emission only, outlined in Table Notice that, in addition to 
the Galactic plane signals, several emissions are noticeable from intensely emitting Galactic 
clouds at intermediate latitudes. Cosmic string searches based on the algorithm presented 
here should be avoided in the highlighted regions. 



ftp://ftp.rssd.esa.int/pub/synchrotron/README.html 
® http : // astrometry.fas .harvard .edu / skymaps 



4. Conclusions 



We studied a direct cosmic string search based on the Cosmic Microwave Background 
(CMB) anisotropies in total intensity, in the context of the existing satehite observations, 
focusing on the case of the ongoing PLANCK survey. The search algorithm uses an alge- 
braic method to detect step like structures in CMB maps, caused by the Kaiser-Stebbins 
effect induced by strings moving on the orthogonal plane with respect to the line of sight. 
We quantified two important aspects towards the application of this algorithm to the ac- 
tual data. 

First, we derived the statistics of null detections due to instrumental noise and purely Gaus- 
sian CMB anisotropies as predicted within the A Cold Dark Matter (ACDM) cosmology. 
We have shown that the biggest cause of confusion for the detection algorithm is caused 
by the CMB acoustic oscillations at wavelengths comparable to the size of the horizon at 
decoupling, while noise and angular resolution which are typical of the operating satellites, 
are less important. The distribution of null detection converges to a x^-distribution. We 
derived the confidence levels of detection of cosmic strings in PLANCK maps by estimating 
the 68%, 95% and 99% of cumulative occurrence of null detection, thresholds give in Table 
III, which corresponds to Gfi ~ 1.5 x 10~^ for 95% of occurrences in terms of cosmic string 
tension, according to the relation between the string tension and temperature shift of CMB 



photons in 2.1 



On the basis of these results, we were able to evaluate the effect of extra-Galactic and 
Galactic foreground emissions on cosmic string searches, by comparison the foreground in- 
duced string signals with the detectability threshold found in the first part of this work. As 
expected, we find that after removal of the brightest sources exceeding 5a times the CMB 
rms, residual confusion due to unresolved point sources does not have a significant impact 
on the string detection algorithm. On the other hand, the diffuse Galactic emission from 
our own Galaxy do cause false detection in a significant fraction of the sky, not because 
of their overall intensity, but rather due to their gradients, reaching intermediate Galactic 
latitudes. By using the current data and models of the diffuse Galactic emissions based 
on off microwave band data as well as the measurements of the WMAP satellite, we deter- 
mine the area of the sky where false detection are expected due to the Galactic emission, 
at the level of 68%, 95% and 99% of cumulative occurrences for each of seven PLANCK 
frequencies, where the thresholds are given in Table j^. We specialize the discussion to seven 
frequency channels on-board PLANCK , constructing for each one the sky mask to be used 
for string searches. The maximum available area in which the Galactic contamination is 
sub-dominant with respect to the contribution from Gaussian CMB fiuctuations is found 
in the frequency interval 70-100 GHz, and represents about 7% of the entire sky. These 
criteria and findings can be used in the actual application of string searches in CMB data. 
For the future work on the pattern search of cosmic strings in CMB temperature data, it 
is desirable to exploit the correlated signals of discrete steps lined up close to each other, 
since the pattern search is mainly pointed at the long strings of cosmological scale rather 
than isolated segments or loops of cosmic strings. The methology provided in this paper on 
finding step signal by a cosmic string segment will be the basis for the long string search. 
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Figure 3: Area overshadowed by strong signals from diffuse foreground emissions for each of the 
7 PLANCK frequencies considered, 30, 44, 70, 100, 143, 217 and 353 GHz. Blue, yellow and red 
colored areas denote the regions where the confusion outcomes by foreground emissions exceed 
68%, 95% and 99% of cumulative occurrences from simulated CMB maps for each corresponding 
frequencies, which are shown in Table 0. 
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